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Abstract 

The micromechanical generalized method of cells model is employed for the predic- 
tion of the effective elastic, piezoelectric, dielectric, pyroelectric and thermal-expansion 
constants of multiphase composites with embedded piezoelectric materials. The pre- 
dicted effective constants are compared with other micromechanical methods available 
in the literature and good agreements are obtained. 


1 Introduction 

Piezoelectric materials are dielectrics which exhibit significant material deformations in re- 
sponse to an applied electric field, and produce electric polarization in response to applied 
mechanical loadings. Piezoelectric materials are used as sensors and actuators, thus form- 
ing important ingredients in smart material systems. Detailed discussion and numerous 
references for piezoelectric sensors and actuators can be found in the book by Gandhi and 
Thompson (1992). A concentric cylinder model that couples elastic, electric, magnetic and 
thermal fields has been presented by Carman et al(1995). 

It is possible to combine piezoelectric materials with passive materials to form various 
types of piezoelectric composite systems. For example one can disperse piezoelectric particles 
into epoxy resin, or embed piezoelectric wires into epoxy resin. Other composite systems 
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Me polymer filled holes in solid piezoelectric ceramic matrices, and porous piezoelectric 
ceramics. Alternatively, one can generate smart materials that consist of fibrous composites 
with embedded piezoelectric actuators and sensors. Discussions of various applications of 
piezoceramic-polymer composites can be found in the book by Moulson and Herbert (1990). 
In all these cases the overall behavior of these composites can be determined by a suitable 
micromechanical model, see Dunn and Taya (1993) and Dunn (1993) for example. Jain and 
Sirkis (1994) considered the issue of damage caused by the growth of voids in piezoelectric 
ceramics. 

The method of cells and its generalization, referred to as the generalized method of cells 
(GMC), is an approximate analytical micromechanical model which is capable of predicting 
the overall behavior of long-fiber, short-fiber and particulate multiphase composites from 
the knowledge of the properties of the individual phases. A review of the method of cells 
followed by a monograph were given by Aboudi (1989,1991). This review has been recently 
updated by Aboudi (1996), where critical assessments of the method and its application by 
various researchers were outlined. As documented, many types of composites (e.g. ther- 
moelastic, viscoelastic, nonlinear elastic and viscoplastic) were analyzed by the method and 
the reliability of the predictions were demonstrated under many circumstances. 

In the present paper, GMC is employed for the prediction of the overall behavior of 
multiphase composites that include one or several types of piezoelectric materials. By im- 
posing mechanical equilibrium and Maxwell’s charge equation in the constituent regions it is 
possible, in conjunction with the continuity of mechanical displacements, tractions, electric 
potential and electric displacements, to establish the electro-mechanical concentration tensor 
that relates the local field in terms of the global one. This readily enables the derivation of all 
effective constants of the multiphase piezoelectric material. The reliability of the predicted 
effective constants is examined by comparisons with other results available in the literature. 

As in previous investigations where GMC was used to develop the macroscale constitu- 
tive equation and supplied as input to analyze composite structures, the derived constitutive 
relations given in this paper can be used to investigate the behavior of composite struc- 
tures with piezoelectric components. Piezoelectric composite plates were recently considered 
by Tauchert(1992) and Jonnalagadda et al(1994). Thus the proposed methodology can be 
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employed to perform thermo-electro-mechanical micro-to-macro analysis of composite struc- 
tures based solely on the properties of the individual phases. 

2 Micromechanical Analysis 

2.1 Model description 

Consider a multiphase composite material in which some or all phases are modeled as thermo- 
piezoelectric materials. It is assumed that the composite possesses a periodic structure 
such that a repeating cell can be defined. In Fig. 1, such a repeating cell is shown which 
consists of N a NpN y rectangular parallelepiped subcells. The volume of each subcell is d Q hply, 
where a, /?, 7 are running indices a = 1, • ■ • , iV Q ; /? = 1, • • • , Np\ 7 = 1, • • ■ , N y in the three 
orthogonal directions, respectively. The volume of the repeating cell is dhl where 

N a N 3 N-, 

d=Y^d a , h = , /=£ *7 (!) 

a=l P=1 7=1 

Any subcell can be filled in general by thermo-piezoelectric materials. Piezoelectric uni- 
directional long-fiber composites, short-fiber composites, porous materials, and laminated 
materials are obtained by a proper selection of the geometrical dimensions of the subcells 
and with an appropriate material fillings. 

The micromechanical model employs a first order expansion of the displacement in the 
subcell (a/3 7 ) in terms of the local coordinates (xf.xf ,xf ) located at the center of the 
subcell. 

Uj M7) = w ia(h ) + *<«> + if Xi (°07) + if V\ (Q/37) » = 1,2,3 (2) 

where are the displacement components at the center of the subcell, and 

•^(“07) and ipi a P^ are the microvariables that characterize the linear dependence of the dis- 
placement on the local coordinates xf ,xf ,xf . In eqn.(2) and the sequel, repeated 

Greek letters do not imply summation. 

Similarly the electric potential (voltage) is expanded linearly in terms of the local 

coordinates of the subcell: 

(^7) = £ 0 (^ 7 > + xf + xf 6 (a/37) + 4 7) 6 M7) (3) 
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( 4 ) 


The components of the small strain tensor are given by 

4° W = + djU i“ W ) i,j = 1,2,3 

where d\ = dj dx[ a ^ , ^ = d/dx^ and d 3 = d/dx^\ 

The components of the electrical field e\°^ in the subcell are obtained from the potential 
f< Q ^> via 

E (*01) = _ d . goto) ( 5 ) 

The constitutive equations that govern the interaction of elastic, electric and thermal 
fields in a piezoelectric medium (that fills the subcell(a/37)) are given by (see Parton and 
Kudryavtesev (1988) for example): 


(a/37) _ M<*Py) (aPy) (apy) p(aPy) \(<*Py) at 
a ij ~ L/ ijkl € kl ~ e kij Z'k n ij LA1 


( 6 ) 


where cr^ 1 ^ are the stress components induced in the material in the subcell, and C \^ ) , 


(opy) 


e «f 7) > and denote the fourth order elastic stiffness tensor, the third order piezoelectric 

tensor and the second order thermal stress tensor of the material in the subcell, respectively. 
The temperature change from a reference temperature is denoted by AT. In addition the 
electric displacement vector is given by the relation 

D ( a py) = e (a/3 7 ) 6 (a/?7) + ^Py) E ^Py) + p Wh) AT (7) 

where and are the second order dielectric tensor and the pyroelectric vector, 

respectively. 

With the linear expansions of the displacements and the potential given by eqns.(2) and 
(3) the static equilibrium of the material within the subcell (afl 7) is satisfied, namely 


A<*P-t) _ 




(8) 


and 


D\f y) = 0 

Let the vectors X^ a/37 ^ and be defined as follows: 


(9) 


jQ(aPy) _ [ Cll , c 22 , e 33 , 2e 23) 2e 13 , 2e 12 , -E lt -E 2l -Erf*" 


( 10 ) 


II I 
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Y<“ w = [ffu.ffa.ffas.ffa.ffu.ffu, (11) 

By excluding the temperature effects in the constitutive equations (6)-(7), one can rewrite 
these equations in the compact matrix form: 

Y^) — z (Q/37) X (a/97) (12) 

where the 9th order coefficient matrix has the following form 


Z( Q / 5 T r ) _ 


(13) 


e ( a Pt) _\( a (h) 

In this equation the square matrix of the 6th-order represents the 4th order stiffness 

tensor, e denotes the transpose of the rectangular 3 by 6 matrix e^ 7 ' that represents 
the corresponding third order piezoelectric tensor, and A^ a/?7) is a square matrix of order 3 
that corresponds to the dielectric tensor. 

The volume average of the stresses and electric displacements D, in the entire repeat- 
ing cell (namely in the composite) is given by 

N a N 0 N-, 


y = E E d a h 0 i y Y<°™ 

aru o=i 0=1 7=1 

and similarly, that of the strains e l3 and electric field components E t are 

1 N a N-, 

^ = 7r,E E E d a h<,l 7 xw 
uru 0=1 0=1 7=1 


(14) 


(15) 


2.2 Interfacial continuity of displacements 

Imposition of displacement continuity conditions at the various subcell interfaces within the 
repeating cell, and at the interfaces between neighboring repeating cells provide the following 
N a Np + N a N-y + NpN y + N a + Np + N y relations between the volume-averaged subcell strains 
and composite strains (Aboudi(1995)): 


Na 

E 

Q = 1 


J (a/37) _ J y 

a Q c n — a en, 


/3=l,---,N0 ; 7= AT 7 


(16) 
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E e 22 ^ = h *22, 
0=1 

a = 1 ,- 

■■,N a 

; 7 = 

1 ,- 

..,at 7 

(17) 

E *7 4f 7 ) =^-33, 

7=1 

a = 1 ,-- 

-,N a ; 

/? = 

1 ,- 

■■,N P 

(18) 

Np N y 

E E W 7) 

0=1 7=1 

= hi C23, 

a = 1 ,-- 

■,n q 

(19) 

No . . . 

E E <w, 4f 7) 

a=l 7=1 

— dl € 13 , 

/? = !,•" 

■,N P 

( 20 ) 

N a Np 

E E ‘If” 

a=l/?=l 

= dh 6 i 2 ) 

7 — 1) ' ' 

-,N . 7 

( 21 ) 


2.3 Interfacial continuity of tractions 

The continuity of tractions at the interfaces (both between subcells and repeating cells) 
provide the following GN a NpNy — (N a Np + N a + NpN y ) — (N a + Np + iV 7 ) conditions. 


(en/3-r) 

a ll 

_ ( Q07 ) 

~ °11 

at= 1,- 

• • , AT Q - 1 ; 0 


,N P ; 

7 = 1 


(22) 

_(®07) 

°22 

_ „(ofa) 
— a 22 

a = 1,- 

■■,n q - 0 = 1, 

.••■.JVJj-i ; 

7 = 1 


(23) 

-(<*0l) 

a 33 

„(<*07) 

~ 033 

a = 1,- 

0 = 1, 

r--,N p 

; 7 = 

I,-" 

,N,-1 

(24) 

-(<*07) 

a 23 

-(<*&7) 

— a 23 

a = 1,- 

• • i N a ] 0 = 1. 

r--,N?-l ; 

7 = 1 


(25) 

(a 

^32 

0n,) = a = 

l,'--,N a \0-- 

= N P ; 

7=1, 

• ’ ■ , N-y — l 

(26) 

_(a^7) 

^13 

-(a07) 

= ff 13 

a = 1,- 

•■,N a - 1; 0 

= 1, • 

,Np ; 

7 = 1 


(27) 
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(a^) _ AaPi) 
°3l — °31 

a = N a ; 0 = 1, - - * , Np ; 

7 = 1 , • ■ • , iV 7 - 1 

(28) 

(a/} 7 ) _ A&0l) 

°\2 — a n a 

= 1,.",/V a -1; 0 = 1, • • 

II 

> 

(29) 

0 — 
— 'C* 

b 

II 

£ 

b 

a = N a ; 0-1 ,Np- 

1 ; 7 = !,-••, N y 

(30) 


where a , 0 and 7 are defined in the following manner, 


a = < 


a + 1 
1 


a < N a 
a = N n 



0 < Np 

0 = Np 


( 7+1 


7 < 


7 = S 


[1 7 = N, 

Note the definition for a , for example, ensures that for a < N Q the neighboring subcell in 
the xx - direction is the one labelled by (a + 1 , /3, 7 ) within the repeating volume element, 
whereas for a = N a the neighboring subcell is within the next repeating volume element 
whose first subcell is (1/37). A similar interpretation for 0 and 7 can be made. 


2.4 Interfacial continuity of electric potential 

The continuity of the electric potential at the interfaces within the repeating cell, and between 
neighboring cells provides N a Np + N a N~, + NpN 7 relations as follows: 


E d ° 

a = l 

= d E u 

0=\,---,Np ; 7 = 1,” 


(31) 


0=1 

= h E 2 , 

r“H 

II 

1— H 

II 

a 


(32) 


7 



£ l 7 = l E 3 , a=l,-,N. ; 0-1, •••,«> (33) 

7=1 

where eqn. (5) and the fact that the average potential gradients are given by 

-2 = &,,= -£ i = 1,2,3 (34) 

OXi 

for all a,P and 7 have been employed in eqns. (31)-(33). 

2.5 Interfacial continuity of the electric displacements 

Finally, the additional 3 N a NpN y - (N a Np + N a N y + N^N-,) required relations are obtained 
by imposing continuity conditions on the electric displacements Di (which ensures that they 
are continuous at the interfaces in the ith directions), that is, 


D{ a0l) 

= d[ 60i) 

Q = 1, 


-l; /? = 

1 

; 7 = 

1, •, Ny 

(35) 

D ( 2 a0l) 

= D[ a0l) 

a-l, 


; /?=!,■■ 

■ , Np 1 

. 7 = 


(36) 

D { ° 0l) 

= Di a ™ 

a = l, 



; 7 

= 1 , • ■ 


(37) 


where a , /3 and 7 are defined as before. 

2.6 Overall constitutive equations 

The combination of the aforementioned continuity relations (16)-(37) forms a system of 
9N a NpN 7 equations in the 9N a NpN^ subcell unknowns. These unknowns can be combined 
and represented by a single vector as follows. 

X s = (X (1U) , X (112) , ..., X^'W) (38) 

This system of equations can be represented in the form 

A X a = K X (39) 
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( 40 ) 


where A is the square matrix of coefficients of order 9N a NpN y , 

X = (en,e22,^33)2e23>2ei3)2ei2,— E\,— E2,—Ez) 

being the average strain and electric field vector, and K is a rectangular matrix of the order 
of 9N a NfiN y by 9, that involves the geometrical properties of the subcells. Note A and K 
consist of the following submatrices: 

Am 0 

A = A g j K = J 

A e J [0 

where A M , involves the elastic mechanical properties of the material, Aq the geometrical 
dimensions of the subcells, Ae the electrical properties, and J the dimensions of the repeating 
cell. 

The solution of eqn.(39) yields the strain and electric field components in all subcells in 
the form 

X s = A X (41) 

where 

A = A 1 K 

and represents the resulting electro-mechanical concentration matrix of the GMC model 

that expresses the local (subcells) variables, X B , in terms of the global (average) variables, 

X. Consequently, eqn.(41) is usually referred to as localization. 

Let the concentration matrix A be partitioned into N a NpN 7 9-order square submatrices 

in the form .. 

A (111) 

A = i (42) 

^{N a N B N^ 

It follows from eqn.(41) that 

X( Q /br) _ X (43) 

Equation (43) expresses the average strain in the subcell (0/^7) in terms of the uniform 
overall strain and electric field (i.e., the applied macrostrain and macroelectric field). 
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Substitution of eqn.(43) into (12) yields 


Y (°07) = 2( q ^A ( q/37) X 


(44) 


Consequently, in conjunction with the averaging procedure given by (14), the following 
effective constitutive relations of the piezoelectric composite can be established 

Y = Z* X (45) 


where the effective elastic stiffness, piezoelectric and dielectric matrix Z* of the composite 
is given by 


N a N 0 N-, 

Z‘ = -7T, £ £ £ <*. h <7 A<«^> 

0/14 a=l 0=1 7=1 

The structure of the square 9 — th order effective matrix Z* is of the form 


(46) 


Z* 


C* e* 4 


(47) 


L e * - A * J 

Thus this micromechanically derived relation readily provides the effective 6 by 6 stiffness 
matrix C*, the effective piezoelectric 3 by 6 matrix e*, and the effective dielectric 3 by 3 
matrix A*. 


2.7 Effective thermal and pyroelectric coefficients 


In establishing the overall behavior of the piezoelectric composite, isothermal conditions were 
assumed (see eqn. (12)). Although it is possible to use the same GMC model to predict the 
effective thermal and pyroelectric constants, it is more convenient to utilize Levin’s (1967) 
result to establish these constants. This approach was previously used by Dunn (1993) to 
establish these constants in conjunction with other micromechanical models. It was also 
utilized by Aboudi (1995) to verify that the effective thermal coefficients and the overall 
inelastic strains as predicted by GMC are in conformity with the corresponding prediction 
by Levin’s method. 

To this end referring to eqns. (6) and (7), let us define the following vector of thermal 
stresses and pyroelectric constants in the material filling the subcell (a0~y): 


j-i(ay37) _ 


A (a/37) A (a/b) \(o/3-y) a (a/37) 
^2 i A 3 1 A 4 ) A 5 


a (a/37) (a/37) (a/37) (a/37) \ 

a 6 >P 1 > P2 iP3 ) 


(48) 



The corresponding global (effective) vector is defined by 


r* 


(AJ.as.asaj.as.as.pJ.pS.pS) 


( 49 ) 


According to Levin’s result, the relation between and T* is given in terms of the 

electro-mechanical concentration matrices AS a ^ as follows. 

1 No. N 0 Ny 

r = [A'“*>] r (so) 

dtl1 a=l 0=1 7=1 

where [A (a/J7) ] r is the transpose of A (q/?7) . The above relation provides the effective thermal 
stress A* and pyroelectric p* vectors of the composite. 

Finally, the effective coefficients of thermal expansion a* and the associated pyroelectric 
constants P* of the composite can be assembled into the vector: 


fl* — (oq, C*2, a 3> a 4> Q 5> a 6> ^1 I ^2 ' ^3 ) 


(51) 


This vector is given by: 


n* = z* -1 r* 


(52) 


It should be mentioned that a similar relation holds for the local quantities, namely for 
the the coefficients of thermal expansion and the associated pyroelectric constants of the 
materials filling the subcells. In such a case Z * and T* should be replaced by Z^ a ^ 7 ^ and 
rC"^), respectively. 


3 Piezothermoelastic Laminated Plate 

The micromechanically established effective constitutive relations 

»,j = c; ju e kl - e- klj E t - A-, AT (53) 

can be utilized to derive the basic equations for piezothermoelastic laminated plates in con- 
junction with the classical plate theory. To this end, let us consider, N layers of orthorhombic 
piezoelectric materials. For this type of material, the constitutive equations, eqns. (53), that 
describes the behavior of a single lamina, can be written in the following form (Tauchert, 
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1992). 


f ' 

On 


C 11 C 12 C 13 0 0 0 


f ' 

l 11 

^22 


c$ 2 C$3 0 0 0 


€22 

<^33 


c 33 0 0 0 


£33 


► — 


4 


^23 


iV 

0 

0 


2^23 

&13 


4 5 0 


2en 

&12 


sym. Cg 6 _ 


, 2e 12 _ 


0 

0 

e 31 




0 

0 

£32 


/ i 





Ei 


0 

0 

eh 

i 

e 2 

► — 4 

0 

e’ 2A 

0 







Ez 


4s 

0 

0 


\ > 


0 

0 

0 



< 


> AT 


(54) 


a; 

A 2 
A3 
0 
0 
0 

When the behavior of the lamina is referred to the plate system of coordinates x,y, z — X3, 
and plane stress conditions are imposed, the following constitutive relations can be obtained 
from (54). 




f ^ 


- 


' 


f ' 



€xx 


0 0 e 3 i 


E x 


Ax 

a vv 

> = [Q) i 

e yy 

> ~ 

0 0 632 

* 

Ey 

— 4 

Ay 

0 X y 


2 e Z y 

k J 


0 0 e 3 g 


E z 

k 


A X y 
k 


> AT 


(55) 


In this equation [0> ; ] are the reduced stiffness coefficients which are given in terms of the 
transformed effective stiffnesses [c^] by the standard relations (see Herakovich (1998)): 


Qij ~ c ij 


c t3 c 3j 

^33 


i,j = 1,2,6 


The transformed piezoelectric constants and thermal stress coefficients are given by 
£31 = £31 c 2 + e^s 2 , £32 = e 31 s 2 + e 32 c , e 36 = (e 31 — e 32 )sc 

A x = AjC 2 4- A 2 s 2 , A y = A*s 2 + AJc 2 , A xy = (AJ - A * 2 )sc 
with c and s denoting the sin and cos of the angle between the x and X\ axis. 
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The stress resultants, N, and moments, M, of the laminate can be readily obtained in 
the form: 

f TVT 1 r A I f fO 1 [ 1M E I I I 

(56) 


N 

M 


> = 


A B 
B D 


6“ 

K 


n e 

M e 


n t 


where A, B and D denote the standard extensional, coupling and bending stiffnesses, and, 
€°and k represent the middle surface strains and curvatures, respectively. In eqn.(56), the 
electric resultants, N E , and moments, M E , are defined by 


(57) 



’ Njf ' 

N 

* \ 
631 

n £ = < 

N y E 

= £ 

^32 | 



fc- 1 



s 

V 


636 1 
* / 



r M E 


r ' 

^31 

m e = < 

X 

M* 

II 

to | *— » * 

632 


K 

J 


636 


> ( E z )k{hl — h\_ j) 


(58) 


where A* denotes the thickness of the kth lamina. The thermal resultants, N T , and moments, 
M T , are given by 


M r = 


Nj 

N 

Af 

-£ 

Nly 

k=l 

Ml ' 

1 " 


4 k= 1 

K , 



Ax 

A. 


A 


xy 


> AT (hk — hfc-i) 


(59) 


i 


v xy 


AT(h 2 fc -fiLi) 


(60) 


4 Results 


In order to assess the reliability of the effective constants of piezoelectric composites predicted 
by GMC, comparisons with other micromechanical methods and with measured data must 
be performed. Dunn and Taya (1993) and Dunn (1993) conducted extensive study of the 
effective constants of piezoelectric composites based on various micromechanical models, and 
performed several comparisons with experimental results. They concluded that the Mori 
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and Tanka (1973) scheme (MT) is the recommended one for application. Consequently, we 
present here comparisons with MT predictions as given by Dunn and Taya (1993) and Dunn 
(1993) in order to verify the validity of GMC predictions. To this end piezoelectric materials 
with hexagonal symmetry in which the axis of symmetry is aligned in the 3-direction are 
considered. The constitutive law of the material that possesses this type of symmetry is 
given by eqn. (54) when 


C 11 — c 22 ’ C 23 — C 13 * c 44 — C 55 > C 66 — 2^ Cl1 C 12 )> e 31 — e 32 > e 24 e 15 ' 


Tables 1-3 present the properties of several piezoelectric materials together with the proper- 
ties of isotropic polymers. 

In Fig. 2 comparisons between GMC and MT predictions are given for the effective 
piezoelectric constants of a porous PZT - 7 A ceramic. Clearly, good agreement is observed. 

Let us define the effective compliance matrix S* to be the inverse of the effective stiffness 
matrix C*. Figure 3 shows the effective elastic compliances 5^ + S* 2 of a piezoelectric 
PZT— 7 A ceramic long fiber reinforcing an epoxy matrix. Here the GMC and MT predictions 
coincide. 

Another set of effective piezoelectric coefficients d £ (which appear when the constitutive 
relation of the piezoelectric material is formulated as strains expressed in terms of stresses 
and electric field) are defined as follows 


In Fig. 4 the above effective piezoelectric constants as predicted by GMC and MT are shown 
for a PZT - 5 ceramic matrix reinforced by long polymer fibers. The agreement is excellent 
for d 33 and satisfactory for d/,. 

Figures 5-7 present the effective coefficients of thermal expansion and the effective py- 
roelectric constants of a polymeric matrix reinforced by long piezoelectric fibers. The pre- 
dictions of GMC and MT coincide, as expected since both micromechanical approaches use 
Levin’s method. 

Figures 8 and 9 present the effective coefficients of thermal expansion and the effective 
pyroelectric constant of a polymeric matrix reinforced by piezoelectric particles. The pre- 
dictions of GMC and MT coincide for the effective thermal properties shown in Fig. 8 and 



closely correspond for the effective pyroelectric coefficient given in Fig. 9. For the same 
composite system, the effective pyroelectric constant —P£/Z£g 1 as predicted by GMC and 
MT is shown in Fig. 10. Here, slight differences can be observed in the vicinity of very high 
volume fraction where steep gradients occur. 

Previous analyses carried out by various researchers (see the update review by Aboudi 
(1996)) have shown that GMC is a reliable and efficient micromechanical model for predicting 
the behavior of various types of composite materials. Based on the present results, it can be 
concluded that GMC is also a reliable micromechanical tool for the prediction of the effective 
parameters of piezoelectric multiphase composites. 

5 Conclusions 

Micromechanical prediction of the effective elastic, piezoelectric, thermal and pyroelec- 
tric constants of piezoelectric composites by the generalized method of cells (GMC) have 
been derived. These constants were verified by comparisons with those predicted using the 
Mori-Tanaka (MT) model and found to be in good agreement. An advantage of the present 
model, over other micromechanical formulations, is its ability to predict the behavior of 
multiphased piezoelectric materials and composites. In this connection it should be em- 
phasized that the Mori-Tanaka method has been previously shown by Ferrari (1991) to yield 
an asymmetric effective stiffness tensor in general, which obviously is a shortcoming of the 
method. Additionally, GMC can be used, to investigate other effects (as previously shown 
by several investigators), such as; the effect of inelasticity of the host material, the effect of 
weak bonding between the constituents, and the effect of fiber distribution and architecture 
on the overall response of piezoelectric composites. Also, the established macromechani- 
cal relations can be readily linked to a structural analysis for the analysis of piezoelectric 
composite structures. 
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Table 1. Elastic material properties. 


Material 

Cn(GPa) 

C n (GPa) 

Cn(GPa) 

C 33 (GP a) 

C 44 (GPa) 

PZT - 7 A 

148 

76.2 

74.2 

131 

25.4 

BaTi0 3 

150 

66 

66 

146 

44 

Epoxy 

8 

4.4 

4.4 

8 

1.8 

Polymer 

3.86 

2.57 

2.57 

3.86 

0.64 


Table 2. Electric material properties. 


Material 

ei 5 {C/m 2 ) 

e 31 (C/m 2 ) 

MC/m 2 ) 

\i(10~ 9 C/Vm) 

A 3 (10 ~ 9 C/Vm) 

PZT - 7 A 

9.2 

-2.1 

9.5 

4.07 

2.07 

BaTi0 3 

11.4 

-4.35 

17.5 

9.86 

11.15 

Epoxy 

0 

0 

0 

0.037 

0.037 

Polymer 

0 

0 

0 

0.079 

0.079 


Table 3. Co efficients of thermal expansion and pyroelectric material properties. 

Material ai(10~ 6 /K) a 3 (10~ 6 /K) P 3 (10 5 N/CK) 

BaTiOz 8.53 1.99 0.133 

Epoxy 60 60 0 
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— x 2 

ng cell in GMC consisting of N a , Np and N Y sub< 
ctions, respectively. 



v f 

Figure 2. — Comparison between the effective piezoelectric 
moduli ejj as predicted by GMC and MT for a porous 
PZT - 7 A ceramic as a function of the porosity volume 
fraction. 
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Figure 3 .— The effective elastic compliances + 

Si*2 of an epoxy matrix reinforced by continuous 
PZT - 7 A fiber as a function of the fiber volume 
ratio. The predictions of GMC and MT coincide. 


Figure 5 . — The effective coefficients of thermal 
expansion a] = «2, <*3 for a continuous BaTi03 
fiber reinforcing a polymer epoxy matrix as a 
function of the fiber volume fraction. The pre- 
dictions of GMC and MT coincide. 


f 




Figure 4. — Comparison between the effective piezoelectric 
moduli d33 and d£ = d3i + d32 + d33 as predicted by 
GMC and MT for a PZT - 5 ceramic matrix reinforced by 
a continuous polymer as a function of the fiber volume 
fraction. 


Figure 6. — The effective pyroelectric coefficients 
P3, of a continuous BaTiC>3 fiber reinforcing a 
polymer epoxy matrix as a function of the fiber 
volume fraction. The predictions of GMC and 
MT coincide. 


ill T 
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Figure 7. — The effective pyroelectric coefficients 
-P^/Z^g -1 of a continuous BaTiC>3 fiber reinforc- 
ing a polymer epoxy matrix as a function of the 
fiber volume fraction. The predictions of GMC 
and MT coincide. 


Figure 9. — Comparison between the effective pyro- 
electric coefficient P3, as predicted by GMC and 
MT for BaTi03 particles reinforcing a polymer epoxy 
matrix as a function of the particles volume fraction. 



Figure 8. — The effective coefficients of thermal 
expansion = <*2, c*3 for BaTi03 particles 
reinforcing a polymer epoxy matrix as a function 
of the particles volume fraction. The predictions 
of GMC and MT coincide. 



Figure 10. — Comparison between the effective 
pyroelectric coefficient -P^/Z^g" 1 , as predicted 
by GMC and MT for BaTiC>3 particles reinforcing 
a polymer epoxy matrix as a function of the 
particles volume fraction. 
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